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Abstract We discuss on noise in heterodyne holography in an off-axis configuration. 
We show that, for a weak signal, the noise is dominated by the shot noise on the 
reference beam. This noise corresponds to an equivalent noise on the signal beam 
of 1 photo electron per pixel, for the whole sequence of images used to build the 
digital hologram. 



1 Introduction 

Demonstrated by Gabor [ 1] in the early 50's, the purpose of holography is to record, 
on a 2D detector, the phase and the amplitude of the radiation field scattered by 
an object under coherent illumination. The photographic film used in conventional 
holography is replaced by a 2D electronic detection in digital holography [2] en- 
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abling quantitative numerical analysis. Digital holography has been waiting for the 
recent development of computer and video technology to be experimentally demon- 
strated [3]. The main advantage of digital holography is that, contrary to holography 
with photographic plates [1], the holograms are recorded by a CCD, and the image 
is digitally reconstructed by a computer, avoiding photographic processing [4]. 

Off-axis holography [5] is the oldest configuration adapted to digital holography 
[6, 3, 7]. In off-axis digital holography, as well as in photographic plate holography, 
the reference or local oscillator (LO) beam is angularly tilted with respect to the 
object observation axis. It is then possible to record, with a single hologram, the 
two quadratures of the object's complex field. However, the object field of view 
is reduced, since one must avoid the overlapping of the image with the conjugate 
image alias [8]. In Phase-shifting digital holography, which has been introduced 
later [9], several images are recorded with different LO beam phases. It is then 
possible to obtain the two quadratures of the object field in an in-line configuration 
even though the conjugate image alias and the true image overlap, because aliases 
can be removed by taking image differences. 

We have developed an alternative phase-shifting digital holography technique, 
called heterodyne holography, that uses a frequency shift of the reference beam to 
continuously shift the phase of the recorded interference pattern [10]. One of the 
advantages of this technique is its ability to provide accurate phase shifts that al- 
low to suppress twin images aliases [11]. This greatly simplifies holographic data 
handling, and improves sensitivity. Moreover, it is possible to perform holographic 
detection at a frequency different from illumination. One can for example detect 
"tagged photons" [12, 13] in ultrasound-modulated optical imaging [14]. One can 
also also image vibrating objects at the frequencies corresponding to vibration side- 
bands [15, 16]. To the end, it is possible to perform Laser Doppler imaging [17] 
within microvessels [18, 19, 20]. 

More generally, our setup can be viewed as a multipixel heterodyne detector that 
is able of recording the complex amplitude of the signal electromagnetic field S 
onto all pixels of the CCD camera in parallel. We get the map of the field over the 
array detector (i.e. £(x,y) where x and y are the pixels coordinates). Since the field 
is measured on all pixels at the same time, the relative phase that is measured for 
different locations (x,y) is meaningful. This means that the field map $(x,y) is a 
hologram that can be used to reconstruct the field S in any location, in particular in 
the object plane. 

In the present paper we will discuss on noise in digital holography, and we will 
try to determine what is the ultimate noise limit both theoretically, and in real time 
holographic experiments. We will see that, in the theoretical ideal case, the lim- 
iting noise is the Shot Noise on the holographic reference beam. In reference to 
heterodyne detection, the reference beam is also called Local Oscillator. We will 
see that the ultimate theoretical limiting noise can be reached in real time holo- 
graphic experiment, by using heterodyne holography [10] in off-axis configuration. 
This combination makes possible to fully filter off the technical noise, whose main 
origin is the LO beam technical noise, opening the way to holography with ultimate 
sensitivity [21, 22]. 
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2 Theoretical noise 

To discuss on noise in digital digital holography, we will consider both the case of 
off-axis holography, where the hologram is obtained from one frame of the CCD 
camera, and the case of phase shifting holography, where the holographic informa- 
tion is extracted from a sequence of M frames. 

We will thus consider a sequence of M frames: 4 to 4f-i (where M = 1 in 
the one shot, off axis case). For each frame 4, let us note h,p,q me CCD camera 
signal on each pixel, where k is the frame index, and p,q the pixel indexes along 
the x and y directions. The CCD signal hp,q is measured in Digital Counts (DC) 
units. In the typical case of the 12 bit digital camera used in experiments below, we 
have < 4,p,o < 4096. For each frame k, the optical signal is integrated by over 
the acquisition time T = 1 // cc d of the CCD camera. The pixel signal h,p,q is thus 
defined by : 




where JJ, p % dxdy represents the integral over the pixel (p, q) area, and where is 
the recording moment of frame k. Introducing the complex representations § and 
<§lo of the fields E and 4lq, we get : 

E{x,y,t) = d?(x,y)e Mr +c.c. (2) 
E W {x,y,t) = S LO {x,y)e^ot +c _ c (3) 

h, P ,q = a 2 ? (K, q \ 2 + KloI 2 + £ P A ■ e {j{ °"- m ^ ] +c.c}) (4) 

where a is the pixel size. To simplify the notations in Eq.4, we have considered that 
the LO field S^q is the same in all locations (x,y), and that signal field S pxl does not 
vary within the pixel (p,q). If S'lo varies with location, one has to replace Sjjo by 
<$LO, P ,q in Eq.4. 

In the single-shot, off-axis holography case, the hologram H is simply H = 4. In 
order to simplify the discussion in the phase shifting digital holography case [9], we 
will consider 4 phases holographic detection (M = 4ri). In that case, the phase shift 
of the LO beam equal to n /2 from one recorded frame to the next. Because of this 
shift, the complex hologram H is obtained by summing the sequence of M frames 
I\ to Im with the appropriate phase coefficient : 

M 

ff=£(;f" 1 4 (5) 

where H is a matrix of pixel H p , CJ , and where M = An in the 4-phases phase-shifting 
case, and M = 1 in the single-shot, off-axis case. We get from Eq.4 : 
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M 

H P ,, = £ U)% P , q = 4na 2 T£ p>q #Z (6) 

k=l 

The complex hologram H p „ is thus proportional to the object field $ p „ with a pro- 
portionality factor that involves S^ Q . 



2.1 The Shot Noise on the CCD pixel signal 

Because of spontaneous emission, laser emission and photodetection are random 
processes, the signal that is obtained on a CCD pixel exhibits a Poisson noise called 
"shot noise". The effect of this Poisson noise on the signal, and on the holographic 
images, is the Ultimate Theoretical Limiting noise, which we will study here. 

We can split the signal h,p,q- we get f° r frame k and pixel (p,q), in a noiseless 
average component (h,p,q) (here { ) is the statistical average operator) and a noise 
component h,p,q'- 

h,p,q = (h,p,q) + h.p.q (7) 

To go further in the discussion, we will use photo electrons Units to measure the 
signal 4 iM . 

We must notice that the local oscillator signal S'lo is large, and corresponds to a 
large number of photo electrons (e). In real life, this assumption is true. For example, 
if we adjust the power of the LO beam to be at the half maximum for the camera 
signal in DC unit (2048 DC for a 12 bits camera), the pixel signal will be about 
10 4 e for the camera used below in experiments, since the "Camera Gain" is 4.8 e 
per DC. This yields two consequences, which simplify the analysis. First, the signal 
h,p,q exhibits a gaussian distribution around its statistical average. Second, both the 
quantization noise of the photo electron signal {h.p.q is an integer in photo electron 
Units), and the quantization noise of the Digital Count signal (Jk,p,q is an integer in 
DC Units) can be neglected. These approximations are valid, since the width of the 
Ik, P ,q gaussian distribution is much larger than one in both photo electron and DC 
Units. In the example given above, (h.p.q) — 10 4 , and this width is ~ 10 2 in photo 
electron Units, and ~ 20 in DC Units. One can thus consider that h,p,q, (h.p.q) an d 
h.p.q are floating numbers (and not integer). Moreover, h.p.q is a random Gaussian 
distribution, with : 

(h.p.q) = (8) 

(i 2 k,,q) = (h,, q ) (9) 

To analyse the shot noise's contribution to the holographic signal H p q , one of 
the most simple method is to perform Monte Carlo simulation from Eq.7, Eq.8 and 
Eq.9. Since h.p.q is ever large in real life (about 10 4 in our experiment), (h.p.q) can 
be replaced by h.p.q (which is measured in experiment) in the right member of Eq.9. 
One has thus: 

(h.p.q) = (h,p,q) - h.p.q (10) 
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2.2 The Object field Equivalent Noise for 1 frame 



90 % 



Fig. 1 1 photon equivalent 
signal (accounting Hetero- 
dyne gain), and shot noise 
on the holographic Local 
Oscillator beam. 



Signal' 

1 photon 
per image 



Detector 
^(CCD Pixel) 



-tof 

Local Oscillator^ 
(LO) 



"shot noise" on LO 
N+/-VN photons 



In order to discuss the effect of the shot noise on the heterodyne signal S p . q S^ 
of Eq.4, let us consider the simple situation sketched on Fig. 1 . A weak object field 
E, with 1 photon or 1 photo electron per pixel and per frame, interferes with a LO 
field Elo with N photons, where N is large (N = 10 4 , in the case of our experiment). 
Since the LO beam signal a 2 T\S'i J o\ 2 is equal to /V photons, and the object field 
signal a 2 T\S p q \ 2 is one photon, we have: 

h,p,q =N+ 1 + i k . p . q + a 2 Tg p . q g£ e ~ +c.c. (11) 

The heterodyne signal S pq S^ is much larger than \S p ^. This is the gain effect, 
associated to the coherent detection of the field $ p „. This gain is commonly called 
"heterodyne gain", and is proportional to the amplitude of the LO field S^ . 

The purpose of the present discussion is to determine the effect of the noise 
term ik,p,q m Eq.ll on the holographic signal H p q . Since H Pt q involves only the 
heterodyne term S p , q S^ (see Eq.6), we have to compare, in Eq.ll, the shot noise 
term ik. P , q , and the heterodyne term S pq S^ . 

Consider first the shot noise term. We have 

(^, 9 } = (W=^+ 1 - i V (12) 

The variance of the shot noise term is thus N l l 2 . Since this noise is mainly related 
to the shot noise on the local oscillator (since Af ^ 1), one can group together, in 
Eq. 11, the LO beam term (i.e. AO with the noise term /j.. p „, and consider that the LO 
beam signal fluctuates, the number of LO beam photons being thus "N ±N l l 2, \ as 
mentioned on Fig. 1 . 

Consider now the the heterodyne beat signal. Since we have N photons on the 
LO beam, and 1 photon on the object beam, we get: 

1 /2 

a 2 T\^lo\ ((« 2 7Kp/) {a 2 T\£ LO \ 2 )) ' =N 1 ' 2 (13) 
The heterodyne beat signal £ VA £t is thus N l/2 = 100. 
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The shot noise term ik,p,q is thus equal to the heterodyne signal S pq S^ corre- 
sponding to 1 photon on the object field. This means that shot noise p q yields an 
equivalent noise of 1 photon per pixel, on the object beam. This result is obtained 
here for 1 frame. We will show that it remains true for a sequence of M frames, 
whatever M = An is. 



2.3 The Object field Equivalent Noise for M = An frames 

Let us introduce the DC component signal D, which is similar to the heterodyne 
signal H given by Eq.5, but without phase factors: 

M 

D=Y,Ik (14) 

k=\ 

The component D can be defined for each pixel (p, q) by : 

M 

D M =£4 )M (15) 

k=l 

Since h,p,q is always large in real life (about 10 4 in our experiment), the shot noise 
term can be neglected in the calculation of D pq by Eq. 15. We have thus: 

M 

Dp.q = £l k , p ,,=Ma 2 T (|# M | 2 + \£lo\ 2 ) (16) 

k=\ 

We are implicitly interested by the low signal situation (i.e. <§ pq <C §lo ) because 
we focus on noise analysis. In that case, the \<§p^ term can be neglected in Eq.16. 
This means that D px/ gives a good approximation for the LO signal. 

M 

D p . q =Y,h,,,^Ma 2 T\£ LO \ 2 (17) 

k=l 

We can get then the signal field l^p.^l 2 from Eq.6 and Eq.17: 

^LL^Ma 2 T\<g p , q \ 2 (18) 

In this equation, the ratio \H Piq \ 2 /D p>q is proportional to the number of frames of the 
sequence (M = An), This means that \H p , q \ 2 /D p iCj represents the signal field \d? P ,q\ 2 
summed over the all frames. 

Let us calculate the effect of the shot noise on \H pXj \ 2 /D p q . To calculate this 
effect, one can make a Monte Carlo simulation as mentioned above, but a simpler 
calculation can be done here. Let us develop \H p q \ in statistical average and noise 
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components (as done for h,p, q in Eq.7): 

H p , q =(H p , q )+h p , q (19) 

with 

An 

h P , q =^j%P,<i (20) 
Let us calculate (\H p ^\ 2 /D p , q ) from Eq.18. Since D pq ~ (D p , q ), we get : 
\H p , q \ 2 \ ^ \{H p , q )\ 2 + (\h p . q \ 2 ) + {{H M )h* M ) + {{H; >q )h M ) 

\A*-) 



D p , q I (D p , q ) 

In Eq.21 the ((H pX! )h* p ) term is zero since h* p is random while (H pxi ) is not. The 
two terms {{H pxi )h* pq ) and ((H* q )h Piq ) can be thus removed. On the other hand, we 
get for \h M \ 2 

An An An 

\ h P,l\ = E +H H ^ Kp,^V,P,q ( 22 ) 

*=1 fc=l k' = \,k!±k 

Since ;'a p c/ and if^, 9 are uncorrelated, the ik, P , q ik> , P , q terms cancel in the calculation 
of the statistical average of |/j p , 9 | 2 . We get then from Eq.9 



An An 

£(I'A VM | 2 ) = £ 
i=l Jt= 1 

Eq.21 becomes thus 



M 2 > - I(l^,?| 2 > = E <W = (Dp, q ) (23) 

l g P,g| 2 \ _ 1(*W1 2 , x (24) 
D M / 

Equation 24 means that the average detected intensity signal (\H p . q \ 2 /D lhq ) is the 

1 /2 

sum of the square of the average object field {\H p , q \) / {(D p ^ q ) 1 ) plus one photo- 
electron. Without illumination of the object, the average object field is zero, and 
the detected signal is 1 photo-electron. The equation establishes thus that the LO 
shot noise yields a signal intensity corresponding exactly to 1 photo-electron (e) per 
pixel. 

The 1 e noise floor, we get here, can be also interpreted as resulting from the 
heterodyne detection of the vacuum field fluctuations [23]. 



3 Reaching the Shot Noise in real life holographic experiment. 



In section 2, we have shown that the theoretical noise on the holographic recon- 
structed intensity images is 1 photo electron per pixel whatever the number of 
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recorded frames is. We will now discuss the ability to reach this limit in real time 
holographic experiment. Since we consider implicitly a very weak object beam sig- 
nal (S <c &lo ), the noises that must be considered are the readout noise of the CCD 
camera, the technical noise from laser amplitude fluctuations on the LO beam, and 
the LO beam shot noise, which yields the theoretical noise limit. 

Consider a typical holographic experiment made with a PCO Pixelfiy 12 bit dig- 
ital camera. The LO beam power is adjusted in order to be at half saturation of the 
digital camera output. Since the camera is 12 bits, and since the camera "gain" is 
4.8e/DC, half saturation corresponds to 2000 DC on the A/D Converter, i.e. about 
10 4 e on the each CCD pixel . The LO shot noise, which is about 100 e, is thus much 
larger than the Pixelfiy Read Noise (20 e), Dark Noise (3 e/sec) and A/D quantiza- 
tion noise (4.8 e, since 1 DC corresponds to 4.8 e). The noise of the camera can be 
neglected, and is not a limiting factor for reaching the noise theoretical limit. 

The LO beam that reaches the camera is essentially flat field (i.e. the field in- 
tensity \<§lo^ lS roughly the same for all the pixels). The LO beam technical noise 
is thus highly correlated in all pixels. This is in particular the case for the noise 
induced by the fluctuations of the main laser intensity, or by the vibrations of the 
mirrors within the LO beam arm. To illustrate this point, we have recorded a se- 
quence of M = An = 4 frames 4 (with k = 0...3) with a LO beam, but without signal 
from the object (i.e. without illumination of the object). We have thus recorded the 
hologram of the "vacuum field". We have calculated then the complex hologram 
H(x,y) by Eq.5, and the reciprocal space hologram H(k x ,k y ) by Fourier transform: 



Fig. 2 Intensity image of 
H(k x ,k y ,0) forM = 4;; = 4 
frames without illumination 
of the object (no signal field 
£). Three kind of noises can 
be identified. Down left (1) : 
shot noise; center (2): techni- 
cal noise of the CCD; left (3) 
: FFT aliasing. By truncating 
the image and keeping only 
the left down part, the shot 
noise limit is reached. The 
image is displayed in arbitrary 
logarithm grey scale. 



H(k Xl k y ) = FFT H(x,y) 



(25) 




The reciprocal space holographic intensity \H\ 2 is displayed on Fig. 2 in arbi- 
trary logarithm grey scale. On most of the reciprocal space (within for example cir- 
cle 1), \H\ 2 corresponds to a random speckle whose average intensity is uniformly 
distributed along k x and k y . One observes nevertheless bright points within circle 
2, which corresponds to (k x ,k y ) ~ (0,0). These points correspond to the technical 
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noise, which is flat field within the CCD plane (x,y), and which corresponds thus to 
low spatial frequency components gathered around the center of the (k x ,k y ) recipro- 
cal space. One see also, on the Fig. 2 image, an horizontal and a vertical bright line, 
which corresponds to k y ~ and k x ~ (zone 3 on Fig. 2). These parasitic bright 
lines are related to Fast Fourier Transform aliases, that are related to the discontinu- 
ity of the signal 4 and H at edge of the calculation grid, in the (x,y) space. 

We have measured (\H\ 2 ) by replacing the statistical average ( ) by a spatial 
average over a region of the conjugate space without technical noise (i.e. over re- 
gion 1). This gives a measurement of (\H\ 2 ), i.e. a measurement of {\H\ 2 ), since the 
space average of \H\ 2 and \H\ 2 are equal, because of the FFT Parceval theorem. We 
have also measured D from the sequence of frames Zj with k — 0...3 (see Eq.14). 
Knowing the camera Analog Digital (A/D) conversion factor (4.8 e/DC), we have 
calculated the noise intensity (\H\ 2 )/(D) in photo-electron units, and we get, within 
10%, one photo electron per pixel for the average noise within region 2, as expected 
theoretically for the shot noise (see Eq.21). 

To verify that we have truly reached the shot noise limit, we have performed a 
control experiment with a camera illuminated by a tungsten lamp powered by a bat- 
tery. The lamp provides here a clean white light source. The lamp voltage is adjusted 
to get half saturation of the camera (about 2000 DC). Like with the laser experi- 
ment described above, we have recorded a sequence of M = An = 4 frames 4 with 
k = 0...3, and we have calculated H(x,y), and H(k x ,k y ). The image of \H(k x ,k y )\ 
we get is very similar to Fig. 2. Moreover, the average noise intensity in region 2 is 
exactly the same as with a laser (one photo electron per pixel). One has thus : 

(\H\ 2 )/(D) = 1 (26) 

This result is expected since the camera "gain" is measured by assuming that the 
noise obtained in clean lamp control experiment is shot noise limited [24]. Assum- 
ing Eq.26, where (\H\ 2 )/(D) depends on the camera "gain", is thus equivalent to 
make a measurement of the "gain" in e/DC Units. The control experiment made 
here, redoes the "gain" calibration made by the camera manufacturer (i.e. PCO). 
We simply get here, within 10%, the same camera gain (4.8e/DC). 



3.1 Experimental validation with an USAF target. 

We have verified that it is possible to perform shot noise limited holography in real 
life, by recording the hologram of an USAF target in transmission. The holographic 
setup is sketched on Fig. 3. We have recorded sequences of M = An = 12 frames, 
and we have reconstructed the image of the USAF target. 

Figure 4 shows the holographic reconstructed images of the USAF target. The 
intensity of the signal illumination is adjusted with neutral density filters. In order 
to filter-off the technical noise, the reconstruction is done by selecting the order 
1 image of the object, within the reciprocal space [8]. Since the 400 x 400 pixels 
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Fig. 3 Setup of the test exper- 
iment with USAF target. L: 
main laser; BS: Beam splitter; 
AOM1 and AOM2: acousto- 
optic modulators; BE: beam 
expander; M: mirror; Al and 
A2: light attenuators. USAF: 
transmission USAF target 
that is imaged. CCD : CCD 
camera. 




CCD 



L ' AOM1 AOM2 



(a) 


(c) 


(e) 


Ei in in i 












(b) 


(d) 


(f) 




' '31 ill t 





Fig. 4 (a,c,d): Reconstructions of an USAF target with different level of illumination 700 (a), 1 (c) 
and 0.15 e/pixel (d). (b): Simulated Shot Noise noise image. (e,f): Simulated reconstructed image 
obtained by mixing image (a) with weight X, and image (b) with weight 1 —X. The weight X is 
1 /700 (e), and 0. 15/700 (f). Images are displayed in arbitrary logarithmic grey scale. 



region that is selected is off-axis, the low spatial frequency noisy region, which 
correspond to the zero order contributions (region 1 on Fig. 2), is filtered-off. 

Figure 4 (a,c,d) shows the reconstructed images obtained for different illumi- 
nation levels of the USAF target. For each image, we have measured the average 
number of photo electrons per pixel corresponding to the object beam, within the 
reciprocal space region that has been selected for the reconstruction (i.e. 400 x 400 
pixels). The images of Fig. 4 correspond to 700 (a), 1 (c), and 0.15 (d) electron per 
pixel for the sequence of M = An = 12 frames respectively (i.e. 700/12, 1/12 and 
. 1 5 / 1 2 e per pixel and per frame) . 

Here, the object beam intensity has been measured by the following way. We have 
first calibrated the response of our camera with an attenuated laser whose power is 
known. We have then measured with the camera, at high level of signal, the intensity 
of the signal beam alone (without LO beam). We have decreased, to the end, the 
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signal beam intensity by using calibrated attenuators in order to reach the low signal 
level of the images of Fig. 2 (a,c,d). In the case of image (a) with 700e/pix, we also 
have measured the averaged signal intensity from the data themselves by calculating 
\H\ 2 /D (see Eq.18). The two measurements gave the same result: 700e per pixel. 

On figure 4 (a), with 700e per pixel, the US AF signal is much larger than the shot 
noise, and the Sinai to Noise Ratio (SNR) is large. On figure 4 (c), with le per pixel, 
the USAF signal roughly equal to the shot noise, and the SNR is about 1 . With 0. 15e 
per pixel, the SNR is low on Fig.4 (d) (about 0. 15), and the USAF is hardly seen. 

It is nevertheless quite difficult to evaluate the SNR of an image. To perform 
a more quantitative analysis of the noise within the images, we have synthesized 
the noisy images of Fig.4 (e,f) by adding noise to the Fig. 4 (a) noiseless image. 
We have first synthesized a pure Noise image, which is displayed on Fig.4 (b). The 
Noise image, which corresponds to the image that is expected without signal, is 
obtained by the following way. From one of the measured frames (for example Iq) 
we have calculated the noise components ik,p,q by Monte Carlo drawing with the 
condition: 

= h,p# (27) 

This condition corresponds to Eq.9 since (Ik,p,q) — h,p,q- We have synthesized the 
image sequence 4 in the following manner: 

h,p,q = h,p,q + h,p,q (28) 

The Shot Noise image of Fig.4 (b) is reconstructed then from the hpq sequence 
with k = 0.. .12 since M = 4n = 12. 

We have synthesized noisy images by summing the noiseless image of Fig.4 
(a) with weight X, with the Noise image of Fig.4 (b) with weight (1 — X). The 
image of Fig. 4 (e) is obtained with X = 1 /700. Figure 4 (e) corresponds thus to 
the same signal, and the same noise than Figure 4 (c) (i.e. le of signal, and le of 
noise respectively). As expected, Fig. 4 (c) and Fig. 4 (e) are visually very similar. 
The image of Fig. 4 is similarly obtained with X = 0.15/700. It corresponds to the 
same Signal and Noise than Figure 4 (d) (i.e. 0.15e of signal, and le of noise). As 
expected, Fig. 4 (d) and Fig. 4 (f) are visually very similar too. 

Here we demonstrated our ability to synthesize a noisy image with a noise that is 
calculated by Monte Carlo from Eq.27 and 28. Moreover, we have verified that the 
noisy image is visually equivalent to the image we have obtained in experiments. 
These results prove that we are able to assess quantitatively the noise, and that the 
noise that is obtained in experiments reaches the theoretical limit of le of noise per 
pixel for the whole sequence of M = An = 12 frames. 



4 Conclusion 

In this paper we have studied the noise limits in off-axis, heterodyne digital holog- 
raphy. We have shown that because of the heterodyne gain of the holographic de- 
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tection, the noise of the CCD camera can be neglected. Moreover by a proper ar- 
rangement of the holographic setup, that combines off-axis geometry with phase 
shifting acquisition of holograms by heterodyne holography, it is possible to reach 
the theoretical shot noise limit. We have studied theoretically this limit, and we have 
shown that it corresponds to 1 photo electron per pixel for the whole sequence of 
frame that is used to reconstruct the holographic image. This paradoxical result is 
related to the heterodyne detection, where the detection bandwidth is inversely pro- 
portional to the measurement time. We have verified all our results experimentally, 
and we have shown that is possible to image objects at very low signal levels. We 
have also shown that is possible to mimic the very weak illumination levels holo- 
grams obtained in experiments by Monte Carlo noise modeling. 
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